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The thermodynamic H — T phase diagram of Bi2Sr2CaCu2 0s was mapped by measuring local 
equilibrium magnetization M(H,T) in presence of vortex 'shaking'. Two equally sharp first-order 
magnetization steps are revealed in a single temperature sweep, manifesting a liquid- solid- liquid 
sequence. In addition, a second-order glass transition line is revealed by a sharp break in the 
equilibrium M(T) slope. The first- and second-order lines intersect at intermediate temperatures, 
suggesting the existence of four phases: Bragg glass and vortex crystal at low fields, glass and liquid 
at higher fields. 

PACS numbers: 74.25. Qt, 74.25. Dw, 74.72. Hs, 64.70. Pf 



The magnetic field vs. temperature (H — T) phase dia- 
gram of the vortex matter in high-temperature supercon- 
ductors, and in Bi2Sr2CaCu208 (BSCCO) in particular, 
has drawn extensive scientific attention The com- 
monly cited thermodynamic phase diagram of BSCCO 
currently consists of a single unified first-order (FO) melt- 
ing line. It separates the low- field quasi long-range or- 
dered Bragg glass (Br GO phase from the high- field liquid 
and glass phases [3 E Gi III El • It is not clear, however, 
whether the two high-field disordered phases are thermo- 
dynamically distinct, or rather reflect a gradual dynamic 
cross-over from liquid into a frozen, pinned state upon 
cooling |3, |8| . In this letter we show that the equilib- 
rium phase diagram of the vortex matter is indeed more 
diverse than the one usually considered. 

Experimentally, one of the main obstacles in mapping 
the low-temperature thermodynamics of the vortex mat- 
ter is its logarithmically slow relaxation rate. Conse- 
quently, the phase diagram has been studied in the past 
mostly through dynamic phenomena. Two prime exam- 
ples are the irreversibility line itself, marking the onset 
of hysteresis, and the second magnetization peak (SMP), 
observed along such hysteretic magnetization loops |9( . 

Recently, though, vortex 'shaking' was shown to be 
extremely effective in catalyzing relaxation at low tem- 
peratures 0, E3 E|. Its application unveiled the in- 
verse melting and the thermodynamic FO transition as 
the phenomenon underlying the non-equilibrium SMP 
0. The 'shaking' method employs the segregated pen- 
etration of an in-plane field component into the highly 
anisotropic BSCCO samples in the form of Josephson 
vortices, which are confined in between the Cu02 planes. 
In the presence of an ac in-plane field, the Josephson 
vortices instantaneously bisect the pancake vortex (PV) 
stacks on their passage, interacting mainly with adjacent 



PVs, while most of the PVs in the stack remain at rest 
0. These occasional interactions agitate pinned PVs, 
assisting them in assuming their equilibrium configura- 
tion. 

Within the present study we performed local magne- 
tization measurements by field and temperature sweeps, 
while utilizing the 'shaking' method, to map the equilib- 
rium phase diagram of the vortex matter. The cross- 
mapping of the FO melting line along both sweeping 
directions shows an excellent agreement. Temperature 
sweeps provided particularly sharp features, with which 
we demonstrate a liquid-solid-liquid sequence of phases. 
Ertas and Nelson 13] have predicted such a liquid-solid- 
liquid sequence to occur within a single temperature 
sweep, but it was never observed experimentally. We 
further find evidence of a novel second order (SO) phase 
transition within the vortex solid phase, which bears im- 
portant consequences regarding the nature of the BrG 
phase. 

The reported results were obtained with a slightly 
over-doped BSCCO crystal with T c « 90 K grown by 
the travelling solvent floating zone method Q This 
specific sample was polished into a triangular prism of 
base 660 x 270 jim 2 and height 70 /mi |l5j (other samples 
yielded similar results, to be presented elsewhere). The 
sample was attached onto an array of eleven 10 x 10 fim 2 
GaAs/AlGaAs Hall sensors. In all measurements taken 
below 60 K the sample was subject to a 10 Hz in-plane 
ac field of amplitude 350 Oe, which was aligned parallel 
to the planes to an accuracy of a few millidegrees. Note 
that according to the anisotropic scaling theory [l| this 
in-plane field is effectively attenuated by a factor 7 ~ 200 
- the anisotropy constant in BSCCO. We found that at 
higher temperatures 'shaking' had no effect on the FO 
transition besides a small broadening (see below). 
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The field sweep mapping of the FO melting line is 
shown in Fig. ^ The collapse of the hysteretic magneti- 
zation into a reversible behavior upon 'shaking' is demon- 
strated in Fig.QJi taken at 32 K. A reversible magnetiza- 
tion step appears instead of the SMP. To better resolve 
the step we plot in Fig.^ the derivative of the measured 
induction with respect to applied field dB/dH. The FO 
transition thus appears reversibly as a 5-like peak (red) 
on top of the dB/dH « 1 (blue) background. Figure 

shows a color scheme of the derivatives dB / dH mea- 
sured by held sweeps within the temperature range 28- 
80 K. The individual melting peaks combine to give the 
locus of the FO transition line H m (T). Its negative slope 
at elevated temperatures becomes positive below 38 K. 
This non-monotonic behavior marks the change in the 
character of the transition from thermally-induced above 
the extremum to disorder-driven in the inverse melting 
region [3, HI3 

In addition to the field sweeps, our experimental setup 
also enables temperature sweeps in the presence of vortex 
'shaking'. It thus allows to measure directly the temper- 
ature dependence of the equilibrium magnetization at a 
constant applied c-axis field. At fields slightly lower than 
390 Oe (e.g. at 380 Oe along the dashed line in Fig. \^a) 
these sweeps should cross the melting line twice. Re- 
markably, the measured local induction in Fig. 03 indeed 
shows two very clear and opposite equilibrium magneti- 
zation steps on both descending and ascending sweeps. 
Note that for clarity we have subtracted from the data a 
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linear slope aT. It is contributed both by the slight tem- 
perature dependence of the Hall coefficients of the sensors 
and by the linear term of the magnitude of the diamag- 
netic equilibrium magnetization, which monotonically 
decreases with temperature. This is the first observation 
of two FO transitions obtained in a single temperature 
sweep. Moreover, the two steps are equally sharp and 
with comparable heights of about 0.15 G. This demon- 
strates that the thermally- and disorder-driven processes, 
responsible for the melting and the inverse melting re- 
spectively, are equivalent mechanisms leading to a FO 
destruction of the quasi-ordered vortex solid. 

The temperatures at which the magnetization steps 
appear along the temperature sweep of Fig. 03 are in 
complete agreement with the melting behavior deduced 
from field sweeps (dotted lines). Therefore, the transition 
line in Fig. [2t is independent of the specific H — T path 
along which it is approached - a mandatory equilibrium 
property. The remaining small hysteresis of 0.1 to 0.2 G 
between downward and upward sweeps apparently results 
from surface barrier effects, while the vortices in the bulk 
are well equilibrated by the 'shaking'. The finite widths 
(about 0.7 K) that the melting steps attain are mainly 
due to a spatial and a temporal averaging mechanisms. 
The first is introduced by the sensor's finite active area 
that averages over the propagating melting front [lf| . 
which results from the spatially inhomogeneous equilib- 
rium magnetization profile [Tjj. The temporal one is a 
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FIG. 1: (color online) (a) 32 K field sweeps with (solid) and 
without (dotted) an in-plane 350 Oe-10 Hz 'shaking' field. 
Reversible steps in magnetization (solid arrows) appear in- 
stead of the hysteretic SMP (open arrows) upon shaking, 
(b) The derivative of the induction with respect to applied 
field dB/dH at 32 K as a color scheme. The first-order tran- 
sition appears as a paramagnetic peak (red) on top of the 
dB/dH « 1 background (blue), (c) Successive mapping of 
the first-order melting line H m (T) measured by field sweeps. 



FIG. 2: (color online) (a) The FO melting line H m (T) mapped 
via field sweeps (open circles), (b) Local induction B(T) in 
the presence of 'shaking', measured upon temperature sweep 
at 380 Oe along the dashed line in (a). A linear slope aT 
was subtracted for clarity. The two equally narrow FO mag- 
netization steps (black segments) show a liquid-solid-liquid 
sequence. The temperatures at which the phase transitions 
occur coincide with those, extracted from field sweeps (doted 
lines). The color code reflects different phases in Fig. 0] 
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by-product of the 'shaking' technique. The in-plane field 
component is known to slightly reduce the melting tem- 
perature 0, 0, ^| . Consequently, our time averaged 
measurement in the presence of the ac in-plane 'shaking' 
field results in an additional broadening due to the in- 
stantaneous periodic shift of the effective local melting 
temperature. 

We thus turn to report the detection of a novel phase 
transition, whose signature is a distinct break in the slope 
of the magnetization M(T). It is visible around 37 K in 
Fig. 03 at 380 Oe, and becomes much more pronounced 
at fields further away from the extremum of the FO melt- 
ing line, as depicted by Fig. 03 The 420 Oe temperature 
sweep (Fig.Et) does not intersect with the FO line, hence 
no steps appear in the local induction. Nevertheless, a 
sharp break in slope is clearly resolved along both de- 
scending and ascending temperature sweeps at T g . A 
sharp reversible break in the induction slope appears also 
in the 350 Oe temperature sweep of Fig. Ob (dotted line) 
in between the two melting steps, hence within the solid 
phase. This non- analytic behavior is emphasized by the 
sharp step in the derivative dB/dT shown in the insets. 
These kinks were found also in other samples and at var- 
ious Hall-sensor locations, and did not depend on the 
sweeping rate. We thus conclude, that this break in slope 
of the equilibrium magnetization M(T) indicates a ther- 
modynamic SO phase transition. 

Mapping of both the first-order H m (T) and second- 
order H g (T) transition lines onto the equilibrium H — T 
phase diagram is given in Fig. 0] The SO line (solid 
dots) intersects the melting curve (open circles), and 
shows weak field dependence throughout the mapped re- 
gion (and therefore cannot be readily observed by field 
sweeps). The resulting phase diagram consists of four 
distinct thermodynamic phases. 
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FIG. 3: (color online) Local induction B(T), measured along 
temperature sweeps while 'shaking'. A reversible sharp break 
in the slope (at the dotted lines) appears both above (a - 
420 Oe) and below (b - 350 Oe) the melting line H m (T) . The 
insets show a corresponding step in the derivative dB/dT, 
signifying a thermodynamic second-order phase transition. 



The high-field part of the novel SO line can be natu- 
rally identified with the long sought glass transition line. 
It asserts that the low-temperature glass phase is indeed 
thermodynamically distinct from the high-temperature 
liquid one. Several experimental studies have observed 
bulk irreversibility features in BSCCO, which appeared 
above H m (T) at about 35 K [H EH- However, all these 
studies probed dynamic or non-equilibrium vortex prop- 
erties. The glass line of Fig. 0] is the first experimental 
evidence of such a thermodynamic transition in BSCCO. 

Yet, the most intriguing result in Fig. 01 is the detec- 
tion of the SO line within the vortex solid region (orange). 
This implies that two distinct thermodynamic phases are 
present in the low-field region below iJ m (T), contrary 
to the common belief that a single BrG phase prevails 
throughout this part of the phase diagram. A number 
of previous studies indicated a depinning line of similar 
topology within the BrG below H m (T) |2£|. However, 
all these measurements probed only the non-equilibrium 
properties, which were consistent with the existing theo- 
retical dynamic predictions and simulations of depinning 
[23. In contrast, the present finding of a thermody- 
namic line requires a more fundamental reconsideration. 

It is interesting to note that in YBa2Cu307 crystals 
a thermodynamic signature of a SO transition within 
the liquid phase has been reported [2^. There, how- 
ever, the SO line emanates from the upper critical point 
of the FO line, directly extending it to higher fields. 
This topology is consistent with several dynamic mea- 
surements in YBa2Cu307 0,HH> alth oug h alternative 
topologies have been also suggested j2f||23. I n contrast, 
our thermodynamic data of BSCCO show that the SO 




FIG. 4: (color online) The thermodynamic phase diagram of 
BSCCO accommodates four distinct phases, separated by a 
first-order melting line H m (T) (open circles), which is inter- 
sected by the second-order glass line H g (T) (solid dots). The 
inset plots an equivalent phase diagram, calculated based on 
Ref. 31, consisting of a second-order replica symmetry break- 
ing lines H g (T) both above (dotted) and below (dashed) the 
first-order transition H m (T) (solid). 
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and the FO transitions are two independent lines that 
intersect each other nearly at a right angle. 

Several theoretical studies have shown that under the 
elastic medium approximation quasi long-range order of 
the vortex lattice is still retained in the presence of 
quenched disorder, giving rise to the BrG phase [j^E^. 
This phase was found to be stable at all temperatures 
(as long as topological excitations are excluded) in sys- 
tems of dimensionality greater than two and lower than 
four. This is probably the reason why a non-topological 
thermodynamic phase transition of the BrG phase was 
hardly ever considered in 3D models. An exception is a 
Josephson-glass line that was suggested to exist within 
the BrG region j^. Still, the general belief is that the 
BrG phase is robust until dislocations proliferate, which 
gives rise to the FO phase transition IHEIII^l- In 2D sys- 
tems, however, the BrG models did find a possible finite- 
temperature depinning transition, above which disorder 
is no longer relevant. Therefore, the observed SO tran- 
sition could be accounted for within these models only 
in the extreme case of vanishing coupling between the 
superconducting layers. 

In contrast, in a recent theoretical work [^J the free en- 
ergy of the vortex matter in the presence of quenched dis- 
order was explicitly calculated under the lowest-Landau- 
level approximation in a 3D model. It was found that two 
transitions are present: a FO melting, at which the quasi 
long-range order is destroyed, and a SO glass transition, 
below which the replica symmetry is spontaneously bro- 
ken. The inset of Fig. ^ shows the phase diagram, calcu- 
lated using a similar effective 2D model with parameters 
optimized for BSCCO [221 . The calculations reproduce 
the measured features very well. Amongst them are the 
melting and inverse- melting behavior, the discontinuity 
of the magnetization slope dM/dT at the glass transition, 
and the H g (T) line itself, which resides at slightly higher 
temperatures as compared to experiment. In addition, 
the calculations show that the portions of the H g (T) line, 
lying above and below the FO transition H m (T) (dotted 
and dashed lines, respectively), are slightly shifted from 
each other. Therefore, they do not cross the melting line 
at a single point, but rather form two closely located 
tricritical points along it. This minute shift can hardly 
be seen in the inset of Fig. 2J and is below our current 
experimental resolution. 

Within this model the high-field glass phase and the 
low-field BrG are strongly pinned and replica symme- 
try broken, whereas the two high-temperature phases are 
replica symmetric and thus reversible. This conclusion of 
reversibility is consistent also with the existing dynamic 
measurements [2(3 . It is therefore tempting to speculate 
that the phase above H g (T) and below H m (T) should 
acquire a true crystalline order. However, since Rosen- 
stein and Li did not calculate the structure factor and our 
measurements do not probe this quantity, the proposed 
vortex crystal phase certainly calls for further experimen- 



tal and theoretical investigations. 

In summary, we present thermodynamic evidence of a 
possibly second-order glass transition line that splits the 
quasi-ordered vortex solid into two distinct phases. By 
comparing the results with existing dynamic measure- 
ments and a new theoretical study, we suggest that the 
two phases are BrG and a vortex crystal. The glass line 
crosses the first-order melting line near its extremum and 
extends to higher fields, giving rise to two thermodynam- 
ically distinct disordered phases - a glass and a liquid. 
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